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Abstract 

The ground-state properties of C20 fullerene clusters are determined in the framework of the Hub- 
bard model by using lattice density-functional theory (LDFT) and scaling approximations to the 
interaction-energy functional. Results are given for the ground-state energy, kinetic and Coulomb 
energies, local magnetic moments, and charge-excitation gap, as a function of the Coulomb re- 
pulsion U/t and for electron or hole doping 5 close half-band filling (\5\ < 1). The role of elec- 
tron correlations is analyzed by comparing the LDFT results with fully unrestricted Hartree-Fock 
(UHF) calculations which take into account possible noncollinear arrangements of the local spin- 
polarizations. The consequences of the spin-density-wave symmetry breaking, often found in UHF, 
and the implications of this study for more complex fullerene structures are discussed. 



PACS numbers: 36.40. Cg, 71.10.Fd, 73.22.-f 



I. INTRODUCTION 



The discovery of the Ceo fullerene^ and the remarkable physical and chemical properties 
resulting from its unique topology and electronic structure have motivated an extraordi- 
nary research activity in past years.— One of the recent questions of interest in this field is 
the possibility of producing even smaller C or Si cagelike clusters and to synthesize novel 
solids using them as building blocks. Thus, several experimental and theoretical stud- 
ies have been performed in order to elucidate the complex mechanisms of formation of 
these nanostructures . 3 i 4 i 5 i 6 i 7 i 8 i 9 In the particular case of C20, which is expected to be the 
smallest fullerene, experiments indicate that the dominant species in laser-vaporization of 
graphite are ring structures, while theory yields different isomers depending on the cal- 
culation method . 1 ^ 1111211 ? 11411 ^ 1 ^ 17 Hartree-Fock (HF) studies predict that the ring is the 
lowest-energy isomer, followed by the bowl (a substructure of the C 60 fullerene) and then 
by the dodecahedral cage . 14 ' 15 ' 16 In contrast, the local density approximation (LDA) to 
density-functional theory (DFT) yields the cage to be more stable than the bowl and the 
ring, whereas in a generalized gradient approximation (GGA) the same ordering as HF is 
obtained . 14 ' 1 ^ 16 These investigations, as well as quantum Monte Carlo (QMC) calculations,— 
indicate that non-trivial electron correlation effects play a central role in the structural and 
electronic properties of these nanoclusters. 

Recently, Prinzbach and coworkers succeeded to produce the cage-structured fullerene 
C20 in the gas phase starting from the perhydrogenated form C2oH2o- 6 This isomer, which 
is not formed spontaneously in carbon condensation or cluster annealing processes, has a 
sufficiently long lifetime. Moreover, photoelectron spectrum (PES) measurements have been 
performed from which its cage-like structure has been inferred.— Theoretical studies by Saito 
and Miyamoto have reproduced the PES thereby confirming the structure assignment given 
in experiment.— In addition, experimental evidence has been provided on the oligomerization 
of C20 fullerenes 8 - and on the formation of a solid phase of C20 dodecahedral which, unlike 
Ceo solids, involve strong C-C bonds between the pentagons of neighboring units. 

Besides its experimental realization, the topology of C 2 o appears to be a particularly 
attractive physical situation for theoretical investigations of correlated itinerant electrons 
in cagelike clusters with pentagonal rings. Similar studies on C12 and Ceo clusters have al- 
ready been performed in relation to alkali-metal-doped C 60 solids like K 3 C 60 and Rb 3 C 60 , 18 
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since one expects that the basic physics behind the superconducting and optical properties 
of these materials should be captured at the scale of an individual molecular constituent. 
The smaller fullerenes have attracted a special attention in this context due to the per- 
spective of achieving larger superconducting transition temperatures than the Cgo-based 
solids. Indeed, a number of theoretical investigations predict a significant enhancement of 
the electron-phonon coupling as the size of the fullerene is reduced from Ceo to C36, C28, and 
finally C2o - 5ll9l2 ° Consequently, theoretical studies aiming at understanding the properties 
of correlated electrons in small fullerene clusters are of considerable interest. 

Lattice models have been used to determine low-energy properties of Ceo which derive 
from the outermost half-filled 7r-electron cloud. In particular, the spin-density distribu- 
tion on the buckyball has been analyzed in a series of papers . 21 ! 22 ! 23 ' 24 ' 25 ' 26 Coffey and 
Trugman 21 calculated the ground-state spin configuration using a classical antiferromag- 
netic (AF) Heisenberg Hamiltonian corresponding to the strongly correlated limit of the 
Hubbard model. They found that the lowest-energy spin structure shows a nontrivial non- 
collinear order which minimizes magnetic frustrations within each pentagonal ring, keeping 
strong AF short-range order in bonds connecting nearest neighbor (NN) pentagons. 21 The 
spin structure in Cqq was also investigated by means of Hubbard or Pariser-Parr-Pople 
Hamiltonians at half-band filling, taking into account on-site and inter-site Coulomb in- 
teractions within the unrestricted Hartree-Fock (UHF) approximation . 23 ' 24 ' 25 ' 26 A common 
result of these investigations is the presence of a magnetic instability for a critical value 
U c of the on-site Coulomb repulsion yielding a magnetic order which resembles that of the 
classical AF Heisenberg model. In addition remarkable noncollinear spin arrangements and 
charge- density redistributions have been obtained in mean-field calculations as a function 
of electron and hole doping. 26 QMC simulations on Ceo at half-band filling and exact diag- 
onalization studies of cagelike C12 support the existence of non-vanishing short-range spin 
correlations . 22 ' 2 ? However, it should be recalled that there is no experimental evidence for 
an spontaneous symmetry breaking. Instead, the spin-density-wave instability should be 
interpreted as an indication of fluctuating spin-spin correlations within the cluster. It is 
therefore interesting to improve on the treatment of correlations in order to quantify the 
possible effects of artificial symmetry breakings on the electronic properties. 

The purpose of this paper is to investigate the ground-state properties of correlated elec- 
trons on a C20 cagelike cluster by using the Hubbard Hamiltonian and a recently developed 
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lattice density- functional theory (LDFT) . 28 ' 29 Previous applications of this approach to ID 
and 2D lattices have yield quite accurate ground-state properties for all band fillings and 
interaction regimes. A particularly interesting feature in the present context is that LDFT 
does not necessarily involve a symmetry breaking in order to account for the effects of elec- 
tron correlation and localization in the strongly interacting limit.— Thus, it appears as an 
appropriate means of improving on UHF calculations. Moreover, besides the aspects specifi- 
cally related to C20 and its topology, the present calculations should be of interest as a finite 
cluster application of new density-functional approaches to lattice fermion models.— 

The body of the paper is organized as follows. In Sec. |H] the model Hamiltonian and the 
main steps in the formulation of LDFT are briefly recalled. Results for the ground-state 
energy, kinetic and Coulomb energies, local magnetic moments, and charge excitation gap 
of C20 are presented and discussed in Sec. IIHI as a function of Coulomb repulsion Uf t and 
for electron or hole dopings 5 close half-band filling (\5\ < 1). Finally, Sec. IIVI summarizes 
the main conclusions and discusses some perspectives of extensions in view of applications 
to lower symmetry structures. 

II. MODEL HAMILTONIAN AND CALCULATION METHOD 

We consider the Hubbard Hamiltonian 31 

H = -t d tcja + U Y1 %™;T> (!) 

on the C20 structure illustrated in Fig. [TJ In the usual notation, c\ a (q ct ) refers to the 
creation (annihilation) operator of an electron with spin a at site i, and n ia = cj^c^ to the 
corresponding number operator. The first term is the kinetic-energy operator with hopping 
integrals = — t < for NNs and = otherwise. The second term takes into account 
intraatomic interactions by means of the on-site Coulomb repulsion integral U. Given the 
lattice structure and the number of atoms N a = 20, the model is characterized by the number 
of electrons N e , or doping 5 = N e — iV a , and by the ratio U/t. 

In order to determine the ground-state properties we consider a density-functional ap- 
proach to electron correlations in a lattice, in which the fundamental variable is the single- 
particle density matrix 7^.— The ground-state energy E gs and density matrix jfj are deter- 
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FIG. 1: Illustration of the geometry of the cagelike C20 cluster. 



mined by minimizing the energy functional 



E[ 1 ]=E K [ 1 } + W[ 1 ] 



(2) 



with respect to 7^-. The first term in Eq. (J2J) is the kinetic energy 



(3) 



ij 

associated with the electronic motion in the lattice. The second term is the interaction- 
energy functional, formally given by22 



where the constrained minimization runs over all iV-particle states ^[7]) that satisfy 



for all i and j. Thus, W represents the minimum value of the interaction energy compatible 
with a given degree of electron derealization or density matrix 7^. The universal functional 
W[7], valid for all lattice structures and hybridizations, can be considerably simplified if the 
hopping integrals are short ranged. For example, if only NN hoppings are considered as in the 
present case, the kinetic energy Ek is independent of the density-matrix elements between 
sites that are not NNs. Therefore, the constrained search in Eq. (|IJ) may be restricted 
to the I ^[7]) that satisfy Eq. ((HJ) only for i = j and for NN ij. Moreover, for highly 
symmetric clusters like C20, one has the same equilibrium value of 7^ for all NN pairs ij, 
and 7jj = n = N e /N a for all sites i. The interaction energy can then be regarded as a 
simple function W(jij) of the density-matrix element between NNs. However, notice that 



W[y] = min U^MlW^n^^}) 



(4) 




(5) 



a 
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restricting the minimization constraints in Eqs. and (0) to NN 7^ also implies that W 
loses its universal character, since the NN map and the resulting dependence of W on 7^ 
are in principle different for different lattice structures . 28 ^ 9 

Two previously investigated scaling approximations to the interaction energy W of the 
Hubbard model shall be used for the calculations . 29 ^ 3 In the first one the functional de- 
pendence is derived from the exact solution of Eq. (0J) for the Hubbard dimer. It is given 



where Eh~f = N a Un 2 /4 refers to the Hartree-Fock energy and gij = (7^ — 7^)/(7^ — 7^) 
measures the degree of electron correlation in a NN bond ij. Here, 7?- > stands for the 
largest possible value of the bond order 7^ for a given N a , band filling n = N e /N a , and 
lattice structure. It represents the maximum degree of electron derealization regardless of 
correlations. 7°° refers to the strongly correlated limit of 7^, i.e., to the largest NN bond or- 
der that can be achieved under the constraint of vanishing W. For half-band filling 7?? = 0, 
while for n 7^ 1, 7°° > 0M Thus, = 1 in an uncorrelated state, and = in the strongly 
correlated limit corresponding, for example, to a fully localized or fully spin-polarized state 
(Nagaoka state). Physically, the scaling hypothesis underlying Eq. (jUJ) means that the rela- 
tive change in W associated with a given change in g 12 is considered as independent of the 
system under study. Exact numerical studies of the functional dependence of W have shown 
that this is a good approximation in a wide variety of ID, 2D, and 3D lattices and band 
fillings. 2 - Notice that £"hf, 7i2, and 7° 2 are system specific. In practice, 7^ is approximated 
by the ferromagnetic fully-polarized 7™ which can be obtained, as -Ehf and 7° 2 , from the 
single-particle electronic structure. 

The dimer approximation to W is very appealing since it combines remarkable simplicity 
and good accurac y 29 > 35 Nevertheless, it also presents some limitations particularly in the 
limit of strong correlations at half-band filling. For N e = N a and small 7^, can be 

expanded as = (l/8)a 2 U'jf j + 0('jf j ) wrt h a 2 = {lij)~ 2 (e.g., a 2 = 2.92 for a pentagonal 
ring and a 2 = 4.16 for C2o)- The exact functional W shows the same behavior but usually 
with a somewhat larger coefficient a ex (e.g., a ex = 3.21 for a pentagonal ring and a 2 = 
5.22 for C 20 )M These discrepancies have direct consequences on the resulting ground-state 
properties for U/t 3> 1, since in this limit 7?? ~ (Az/a)(t/U) and E gs ~ —(2z 2 /a)(t 2 /U), 
where z is the local coordination number (Ek = —zt / ~f i j). In order to improve the accuracy a 
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more flexible approximation has been proposed which is also based on the scaling properties 
of W, and which recovers the exact dependence on 7^' for 7^ — > 0. In this case the interaction 
energy includes a fourth-order term in and is given by 



where k = a ex /a2- In a pentagonal ring n = 1.10, while in the C20 cluster k = 1.26. 
These values should be compared with the dimer result re = 1, for which Eq. (JJJ) reduces to 
Eq. (jnj). Thus, the 4th-order term appears as a moderate correction to the dimer or 2nd- 
order approximation (gfj < 1). Previous applications to ID and 2D systems have shown that 
Eq. ((7j) provides a systematic improvement on the ground-state properties for all [//£.— In 
the following section, Eqs. © and (0) are applied to determine several electronic properties 
of the C20 Hubbard cluster in the framework of LDFT. 

III. RESULTS 

The model is characterized by the dimensionless parameter U/t and by the doping 5 = 
N e —N a . We consider the complete range of repulsive interactions (U > 0) and single electron 
or hole dopings (\5\ < 1). Commonly accepted values of U/t for the 7r-electrons in carbon 
fullerenes correspond to the intermediate range 2 < U/t < 5 . 23 i 24 i 25 i 26 Before discussing the 
results for the C20 cluster it is useful to consider a single pentagonal ring as a preliminary 
test on the accuracy of the method, since the pentagon constitutes the basic building block 
of the C20 structure, and since it is small enough so that exact diagonalizations can be 
easily performed. Fig. 121 shows results for the ground-state energy E gs of the N a = 5 ring 
as a function of U/t. In the uncorrelated limit E gs decreases monotonously with increasing 
number of electrons N e — N a + 5 < 6 since the three lowest single-particle eigenvalues of 
the ring are negative. On the other side, for large U/t, E gs vanishes at half-band filling 
while for negative (positive) doping E gs (E gs — U5) remains negative. In fact, for 5 = the 
suppression of charge fluctuations always implies electron localization, whereas for \5\ — 1 a 
finite kinetic energy persists, even for U/t ^> 1, due to the derealization of the extra electron 
or hole. The LDFT calculations reproduce remarkably well the crossover between weak and 
strong correlations as given by the exact diagonalizations. Appreciable discrepancies are 
found only for 5 = — 1 and U/t > 6, in which case the largest relative error amounts to 




(7) 
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U/(U+4t) 

FIG. 2: Ground-state energy E gs of the Hubbard model on a pentagonal ring as a function of 
Coulomb interaction U/t for different numbers of electrons N e = N a + 5 close to half-band filling 
(N a = 5). For 5=1, results for E gs — U are shown. The solid curves correspond to lattice density- 
functional theory (LDFT), the crosses to exact numerical diagonalizations, and the dashed curve for 
5 = to unrestricted Hartree-Fock (UHF) calculations including noncollinear spin arrangements. 

\ E g S DFT ~ E gs\/\ E gs\ - °- 07 for & 7* - 12 - In the other cases ( S = or x ) the LDFT results 
are practically indistinguishable from the exact ones. 

In order to quantify the effects of electron correlations we have also performed calculations 
by using the fully unrestricted Hartree-Fock (UHF) approximation which corresponds to the 
most general single-Slater-determinant wave-function and which allows for noncollinear site- 
dependent spin polarizations (Si).— Figure 121 shows the UHF results for E gs of the pentagonal 
ring at half-band filling (5 = 0). For U/t < U c ~ 0.27 the only solution to the self-consistent 
equations is non- magnetic (i.e., (Si) = for all i), while for U/t > U c local magnetic mo- 
ments (Si) set in. These increase monotonously with U/t reaching saturation in the limit of 
U/t — » oo. Concerning the magnetic order one observes that the (Si) are all coplanar and 
that the angle between NN moments is 47r/5 for all U/t > U c . This amounts to split one 
parallel-spin frustration among five bonds which corresponds, as expected for half-band fill- 
ing, to the solution of a classical Heisenberg model with antiferromagnetic NN interactions. 
One observes that UHF reproduces the U/t dependence of E gs qualitatively well, including 
the strongly correlated limit where E gs — ► 0. This is achieved by localizing the electrons 
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FIG. 3: (a) Kinetic energy Ek, (b) Coulomb energy Ec, and (c) local magnetic moments (Sf) 
of the Hubbard model on a pentagonal ring at half-band filling (N e = N a = 5). The dashed- 
dotted and full curves correspond to lattice density-functional theory (LDFT) using, respectively, 
Eqs. (jHJ) and Q for the interaction-energy functional. The dashed curves refer to the unrestricted 
Hartree-Fock (UHF) calculations including noncollinear spin arrangements, and the crosses to exact 
numerical diagonalizations. 

through the formation of increasingly large local magnetic moments as U/ 1 increases. Never- 
theless, significant quantitative discrepancies with the exact numerical solution are observed 
already for U/t > 2, which reflect the importance of electron-correlation effects. These can 
be regarded in part as the result of fluctuations between different spin configurations which 
cannot be recovered in a single-determinant state.— In contrast, LDFT gives a much bet- 
ter description of the ground-state, particularly in the intermediate and strong interacting 
regimes. UHF overestimates E gs appreciably, while the differences between exact and LDFT 
results are very small. 

More detailed information on the correlation effects and on the accuracy of the calcula- 
tions is obtained from the kinetic energy Ek, Coulomb energy Ec, and local square magnetic 
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moments (Sf) presented in Fig. El LDFT accurately reproduces the exact numerical solution, 
showing that the very good results obtained for E gs are not the consequence of important 
compensations of errors. Moreover, the differences between 2nd-order and 4th-order approx- 
imations to W [see Eqs. (JHJ) and (J7J)] are quite small, typically of the order of the differences 
between the 4th-order calculations and the exact results. Nevertheless, one may observe 
that the interaction-energy functional provides a consistent improvement with respect 
to W^. In fact taking into account the 4th-order term yields a small increase (reduction) 
of Ec for U/t < 4 (U/t > 4) as well as a small reduction of \Ek\ for U/t > 4 which bring 
the outcome of LDFT systematically closer to the exact solution. 38 

In contrast to LDFT the UHF calculations do not provide a correct description of the 
kinetic and Coulomb energies separately. On the one side, for U/t < 4, UHF overestimates 
the Coulomb energy Ec due to an underestimation of the local moments (Sf) (see Fig. OJ). 
On the other side, for U/t > 4, both Ec and \Ek\ are appreciably underestimated. The 
inaccuracies in Ec and \E%\ for U/t > 4 tend to cancel each other, which improves to some 
extent the UHF result for E gs (see Fig. EJ). Notice that UHF reproduces quite accurately 
the local moments for U/t > 4 (see Fig. EJ). However, a single-determinant state fails 
to incorporate correlations, and in particular the fluctuations between different equivalent 
spin configurations, which should restore the broken spin symmetry and which are often 
important in finite systems.— This reflects the limitations of the broken symmetry UHF 
solutions, even in the most general noncollinear case.— Notice that in LDFT no artificial 
symmetry breaking is required in order to account for correlation-induced localizations and 
the resulting effects on the ground-state properties. The present LDFT formalism with 
the dimer or 4th-order scaling approximations to interaction-energy functional describes 
correctly the electronic correlations in a pentagonal Hubbard ring. It can be therefore 
expected that the method should remain accurate in more complex systems with five-fold 
symmetry such as the C20 cluster to be discussed below. Moreover, comparison with exact 
results also shows that the differences between UHF and LDFT results for E gs , Ek, and Ec 
are a good quantitative estimation of the correlation energies and their U/t dependence. 

In Figs. 13] and 03 LDFT and UHF results are given for several ground-state properties 
of the Hubbard model on the cagelike C 20 cluster illustrated in Fig. [TJ In the uncorrelated 
limit the ground-state energies E gs are the same for \8\ < 1, a consequence of the degeneracy 
of the single-particle spectrum. As U/t increases, E gs increases monotonously with a slope 
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U/(U+4t) 

FIG. 4: Ground-state energy E gs of the Hubbard model on the cagelike C20 cluster illustrated 
in Fig. ^ Results are given as a function of Coulomb interaction U ft for different numbers of 
electrons iV e = N a + 5 close to half-band filling (N a = 20). In the inset E gs — U is shown for 
5 = 1. The solid curves correspond to lattice density- functional theory (LDFT) using Eq. (JJJ) for 
the interaction-energy functional. The dashed curves refer to unrestricted Hartree-Fock (UHF) 
calculations including noncollinear spin arrangements. 

dEg S /dU = J2i( n ii n ii) > 0. Since the average number of double occupations is larger for 
larger number of electrons the degeneracy with respect to 5 is removed at finite U/t. On 
the other side, in the strongly correlated limit, E gs vanishes at half-band filling, while for 
5 = — 1 (5 = 1) E gs (E gs — U) remains negative. The U/t dependence of E gs in C 2 o follows 
similar trends as in the pentagonal ring. The UHF calculations reproduce the exact limits for 
U/t = and U/t = 00, and agree qualitatively with LDFT for finite values of the Coulomb 
repulsion. However, the quantitative differences can be quite significant particularly for 
intermediate and large U/t (e.g., 26% for U/t ~ 36). 

The effects of correlations are more significant in the kinetic energy Ex, Coulomb energy 
Eci and local square moment (Sf) shown in Fig. El for 5 = 0. The differences between 
UHF and LDFT calculations for Ek and Eq are similar to the corresponding differences 
between UHF and exact results for the pentagonal ring and can be qualitatively understood 
in similar terms. In the case of C20 one may notice that the UHF results for the focal square 
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FIG. 5: (a) Kinetic energy Ek, (b) Coulomb energy Ec, and (c) local magnetic moments (Sf) of 
the cagelike C20 Hubbard cluster at half-band filling (6 = 0). The dashed-dotted and full curves 
correspond to lattice density-functional theory (LDFT) using, respectively, Eqs. © and ((7J) for 
the interaction-energy functional. The dashed curves refer to unrestricted Hartree-Fock (UHF) 
calculations including noncollinear spin arrangements. 

moment (Sf) are the same as in the uncorrelated limit for U/t < 2.2. Therefore, they are 
significantly underestimated. For stronger Coulomb repulsions (U/t > 2.2) (Sf) increases 
rapidly reaching values that are very close to the LDFT results for U/t ~ 4 (see Fig. |3J). 
This behavior is related to the U/t dependence of the self-consistent solution of the UHF 
equations. For < U/t < 3.1 the UHF charge distribution is not uniform due to the four- 
fold degeneracy of the single-particle spectrum at the Fermi energy. Moreover, for U/t < 2.2 
the UHF ground state shows a weakly ferromagnetic solution with very small local spin 
polarizations. In this way advantage is taken from the single-particle degeneracy in order 
to reduce the Coulomb-repulsion energy without increasing significantly the kinetic energy. 
For U/t > 2.2 the system starts to develop an AF spin-density wave which coexists with a 
weak charge- density wave for 2.2 < U/t < 3.1. Only for U/t > 3.1 the AF-like magnetic 
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FIG. 6: Illustration of the distribution of the local spin polarizations (Si) in an undoped cagelike 
C20 cluster as obtained by using the Hubbard model and the unrestricted Hartree-Fock approxima- 
tion ((5 = 0). Due to spin-rotational invariance only the relative orientations of (Si) are physically 
meaningful. The three-dimensional cluster structure (see Fig. ^) is mapped onto a plane in order 
to ease the visualization. The radius of the circle on each atomic site i is proportional to |(Si)|. 
The arrows represent the projection of (Si) onto the xy plane which is the plane containing the 
outermost pentagon. Shaded circles (open circles) indicate that the perpendicular component (Sf) 
is positive (negative). 

order valid for strong interactions is definitely settled and the charge distribution becomes 
then uniform throughout the cluster. An illustration of the UHF magnetic order in the 
large U / 1 regime is shown in Fig. |UJ Notice the noncollinearity of the local spin polarizations 
(Si), which form angles of about 138 degrees between NNs. This value is smaller than the 
angle between NN (Si) in an isolated pentagonal ring (144 degrees) because of frustration 
effects between the pentagons in the C20 topology (see Figs. [T] and EJ). Further increase of 
the Coulomb repulsion beyond U/t — 3.1 only results in an increase of the size of the local 
moments keeping the magnetic order unchanged. 
The charge excitation gap is given by 

AE C = E gs (N e + 1) + E gs (N e - 1) - 2E gs (N e ) , (8) 

where E gs refers to ground-state energy. It is a property of considerable interest since it 
measures the low-energy excitations associated with changes in the number of electron N e 
and is thus very sensitive to electron correlation effects. In Fig. [7| results are shown for AE C 
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FIG. 7: Charge excitation gap AE C of the Hubbard model on the cagelike C20 cluster as a 
function of Coulomb interaction U/t. The solid (dashed) curves correspond to LDFT (UHF). In 
the inset figure results are given for the pentagonal Hubbard ring. Here the crosses refer to exact 
diagonalizations. 

of the C20 Hubbard cluster at half-band filling (N e = N a ). The corresponding calculations 
for the pentagonal ring are reported in the inset for the sake of comparison. At half-band 
filling AE C = for U = 0, despite the finite size, due to the degeneracy of the single-particle 
spectrum in both the ring and C 2 o- AE C increases with increasing U/t approaching the 
limit AE —>([/ — tu 6 ) for U/t — > 00, where Wb is the bottom of single-particle spectrum of 
the cluster = 2t for the ring and Wb = St for C 2 o)- The LDFT results and the exact 
calculations for N a = 5 yield a linear dependence of AE C as a function U/t for U/t <C 1. This 
implies, in particular for C20, that AE C > for arbitrary small U/t, in contrast to the UHF 
result. Comparison with exact diagonalizations for N a = 5 indicates that the LDFT gap is 
not overestimated. The UHF calculations fail to reproduce the correct U/t dependence of 
AE C showing once more the importance of correlations. Notice that the differences between 
UHF and exact or LDFT results are qualitatively similar for C 2 o and for the pentagonal 
ring, which indicates that the contribution of correlations to AE C are comparable in both 
cases. 
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IV. CONCLUSION 



A density-functional approach to lattice-fermion models has been applied to study the 
ground-state properties of the Hubbard Hamiltonian on a cagelike C20 cluster. The ground- 
state energy, the kinetic and Coulomb contributions, and the charge excitation gap have 
been determined as a function of U/t for dopings 5 close to half band filling. The impor- 
tance of electron correlations has been quantified by performing noncollinear Hartree-Fock 
calculations. Comparisons with exact diagonalizations and with UHF in the case of a pen- 
tagonal ring demonstrate the ability of LDFT to describe subtle correlation effects in a very 
simple and systematic way. It is therefore expected that LDFT should be an efficient tool 
for investigating the properties of strongly correlated fermions in more complex clusters and 
nanostructures, which are otherwise very difficult to tackle with other accurate methods. 

The present LDFT study of clusters encourages the development of more flexible 
interaction-energy functionals in order to determine the properties of finite systems with 
lower symmetry. As a first step it would be interesting to extend the calculations to fullerene 
structures having different types of bonds, for example Cgo- In these cases the functional 
dependence of W[7y] is more complex and one has to consider at least two different NN 
density-matrix elements, within and between pentagons, even if one assumes that the hop- 
ping integrals are the same. The symmetry of these bonds is actually different despite the 
fact that all sites are equivalent. The interaction-energy functional should be generalized 
following previous works on dimerized ID chains.— The domain of representability of 7^ and 
the accuracy of the scaling hypothesis could be explicitly tested for smaller similar clusters 
like C12 by solving Levy's constrained search . 28 i 32 Furthermore, the functional dependence 
of W for low-symmetry clusters deserves to be explored by taking into account charge trans- 
fers and eventual redistributions of the spin density between non-equivalent sites. This 
should open the way to investigations of correlation effects in a variety of complex physical 
situations such as metal clusters, disordered systems, and low-dimensional nanostructures. 
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